Indication on the process-dependence of the Sivers effect 
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We analyze the spin asymmetry for single inclusive jet production in proton-proton collisions 
collected by AnDY experiment and the Sivers asymmetry data from semi-inclusive deep inelastic 
scattering experiments. In particular, we consider the role color gauge invariance plays in determin- 
ing the process-dependence of the Sivers effect. We find that after carefully taking into account the 
initial-state and final-state interactions between the active parton and the remnant of the polarized 
hadron, the calculated jet spin asymmetry based on the Sivers functions extracted from HERMES 
and COMPASS experiments is consistent with the AnDY experimental data. This provides a first 
indication for the process-dependence of the Sivers effect in different processes. We also make predic- 
tions for both direct photon and Drell-Yan spin asymmetry, to further test the process-dependence 
of the Sivers effect in future experiments. 
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The investigation of nucleon's sub-structure has en- 
tered a new era. In past decades an understanding of 
nucleons in terms of quarks and gluons (partons), the de- 
grees of freedom of Quantum Chromodynamics (QCD), 
has been successfully established. Progress was achieved 
in constructing a "one-dimensional" light-cone picture of 
the nucleon based on the longitudinal motion of par- 
tons in fast moving nucleons. In recent years theoreti- 
cal breakthroughs extended this description in the trans- 
verse as well as light-cone momentum space (three di- 
mensions). Transverse-spin dependent observables, such 
as single transverse spin asymmetries (SSAs) provide firm 
evidence for a three-dimensional tomography of the nu- 
cleon due to a non-trivial correlation between the trans- 
verse spin and the parton's transverse momentum, and 
present unique opportunities to study QCD dynamics, 
particularly QCD factorization and universality of the 
parton distributions 

Large SSAs have been measured in fixed-target and 
collider mode in single inclusive particle production in 
nucleon-nucleon scattering experiments [2j and semi- 
inclusive deep inelastic lepton-nucleon scattering (SIDIS) 
experiments [3j, |4|. Two different yet related QCD fac- 
torization formalisms have been proposed to describe the 
asymmetries. One relies on the so-called transverse mo- 
mentum dependent (TMD) factorization [f^ @| , which is 
valid for the processes with two characteristic scales; for 
example the photon's virtuality Q and Ph± of the pro- 
duced hadron in SIDIS, where Aq CD < P 2 ± < Q 2 . In 
this formalism transverse spin effects are associated with 
naive time-reversal-odd TMD parton distribution func- 
tions and fragmentation functions (PDFs and FFs). On 
the other hand is the collinear factorization formalism at 
next-to- leading power (twist-3) in the hard scale jzLUJ. 
This approach is valid for processes with only one char- 
acteristic hard scale, for instance, the transverse momen- 
tum Pfi ± 3> Aq CD of the produced hadron in proton- 



proton (pp) collisions. It describes the spin asymme- 
try in terms of twist-3 three-parton correlation functions. 
One of the well-known examples is the so-called Efremov- 
Teryaev-Qiu-Sterman (ETQS) function T q ,p(x,x) Q- 

Of central importance in the study of SSAs is the 
Sivers jl[ effect which has attracted great attention in 
recent years. In part this is due to the unique prediction 
from TMD factorization theorems that the Sivers effect 
is process-dependent: that is its existence relies on the 
initial-state and final-state interactions (ISIs and FSIs) 
between the struck parton and the remnant of the polar- 
ized hadron. These interactions depend on the color flow 
of the specific scattering process considered, thus giv- 
ing rise to process-dependent Wilson lines in the gauge- 
invariant definition of the relevant TMD PDFs - in this 
case so-called Sivers functions fyp{x, k\). The often dis- 
cussed case is the difference between the FSIs in SIDIS 
and the ISIs in Drell-Yan (DY) production in pp colli- 
sions which leads to an opposite sign in the Sivers func- 
tion probed in these two processes, indicating that the 
Sivers function is not universal flfjj ]. 

On the other hand in the twist-3 collinear factoriza- 
tion approach, the process-dependence of the ISIs and 
FSIs is absorbed into the short-distance, perturbative cut 
scattering amplitudes, where the relevant twist-3 three- 
parton correlation functions are universal. As a result, 
TMD and collinear twist-3 factorization formalisms are 



closely related to each other The relevant functions 
- the Sivers function and the ETQS twist-3 function are 
connected through the following relation 12, 13| . 



T Qi f(x, x) = — / d k 



^/^(x^DlsiDIS- (1) 



where the subscript emphasizes that the Sivers function 
is probed in the SIDIS process. In other words, starting 
from the Sivers functions extracted from SIDIS, one can 
derive a functional form for ETQS function T qy p(x, x). In 
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combination with the calculable short-distance cut scat- 
tering amplitudes, one should be able to predict the SSAs 
of inclusive particle production in pp collisions. How- 
ever, a recent study [13l4l5| for inclusive hadron produc- 
tion in pp collisions shows that such calculated SSAs are 
opposite to those measured in the experiments. This is 
known as the "sign mismatch" problem. Whether this 
finding reflects the inconsistency of our theoretical for- 
malism is a very important question and needs to be ex- 
plored both theoretically and experimentally. However, 
since the SSAs of inclusive hadron production can also re- 
ceive contributions from the fragmentation process [t| , a 
thorough analysis demands including such contributions. 

A new opportunity presents itself however, with a re- 
cent inclusive jet measurements performed at the AnDY 
experiment at RHIC [l6j . Since the jet spin asymmetry 
does not involve fragmentation contributions, this paves 
the way to precisely test the process-dependence of the 
Sivers effect in different processes as well as explore the 
consistency of the TMD and collinear twist-3 factoriza- 
tion formalisms. This is the main purpose of our paper. 
We analyze the spin asymmetry for single inclusive jet 
production in pp collisions collected by the AnDY exper- 
iment and the Sivers asymmetry data from SIDIS experi- 
ments. We assess whether they are compatible with each 
other; in other words, whether the jet asymmetry is con- 
sistent with our expectation on the process-dependence 
of the Sivers effect. 

We start with the basic formalism for the SIDIS SSA. 
For hadron production in SIDIS at low transverse mo- 
mentum P hX , e(£) + A^(P,s x ) -> e(l') + h{P h ) + X, 
within the TMD factorization formalism, the differential 
cross section for the Sivers effect reads [171 ]. 

da 



Fuu,t + sin( 



)F, 



sm(tj> h -ip a ) 



UT.T 



(2) 



dV.S. 

where phase space dV.S. = dxBdyd4> s dzhd(f>hPh^dPhi^ 
with the standard SIDIS kinematic variable xb, y, and 
z/j. The normalization factor cr = ^(^b, V, Q 2 ) an d 
the structure functions Fuu.t and ~~ are de- 

fined in Ref. [14| . The Sivers asymmetry measured in 
the experiments is defined by A™,^ h ~'^ (x B ,z h ,p h± ) = 

&0\XB,y,Q )i t u TT /0-o(x B ,y,Q )riJU,T- 

On the other hand, the single inclusive jet produc- 
tion in transversely-polarized pp collisions, A(Pa,s±) + 
B(Pg) — > jet(Pj) + X, only receives the Sivers type 
of contributions. Within the collinear factorization for- 
malism, the spin-dependent differential cross section 
dAa(sj_) = [da(s±) — da(— s±)] /2 can be written as, 



E 



dAa(sj 



a 2 ^ 



a . b 



dx dx' 

——rfb/B{x ) 



T a ,F{x, x) - x—T a F (x, x) 
dx 



x^H!™ s (s,i,u)S(s + i+u), (3) 



where b runs over all parton flavors, fb/B( x ') * s the 
collinear PDF in the unpolarized proton, and s, t, and u 
are the standard partonic Mandclstam variables [14|, [l8j . 
Hab^>c represent the cut scattering amplitudes for the 
partonic process ah — > cd with the expressions given 



15l . Il8j . They incorporate both the ISIs and FSIs 



between the active parton and the remnant of the po- 
larized proton; thus the process-dependence of the Sivers 
effect for jet production is included in -fff ^™^ s . The SSA, 
An is computed from the ratio of the spin-dependent to 
the spin-averaged cross section, 



A N = E. 



dAo~(sj 



Ej 



da 
d?P~j 



(4) 



where the spin- averaged differential cross section Ej °" 
in the denominator is defined in Ref. 1141. 
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FIG. 1. Description of the HERMES 0] data for tt+ pro- 
duction as a function of Bjorken xb- The solid lines are the 
central values from Tableland the shaded region corresponds 
to the parameter scan as explained in the text. 

whether the inclusive jet data in pp collisions are consis- 
tent with the Sivers asymmetry data in SIDIS processes, 
we perform a global fit of the SIDIS Sivers asymmetry 
data collected by the HERMES and COMPASS exper- 
iments 0, |3] to extract the Sivers functions. We then 
derive the functional form for twist-3 ETQS function 
Tq t p(x,x) with the help of Eq. ([TJ and in turn compute 
the jet spin asymmetry An from Eq. ((4]) to be compared 
with the data collected by AnDY experiment jl6| . 

We adopt the Gaussian forms in Ref. [l!| for the spin- 
averaged PDFs, /a/A^^l) and FFs D h / a (z 7 p^), with 
the Gaussian width, (k ± ) = 0.25 GeV 2 and = 0.2 
GeV 2 . The quark Sivers function f^(x,k 2 ± ) for SIDIS 
is parameterized as, 



(5) 



where the fc^-dependence h(k±) = Vtejfe k ±/ M i J 
with M the proton mass, and the ai-dependent coefficient 
N q {x) = N q x a "{\- xf"(a q + fi q ) aq+fiq /(q,"'^/')- For 
the purpose of our fit, we use GRV98LO for the spin- 
averaged collinear PDFs 2l| and DSS parametrization 
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for collinear FFs [22J. During the fit we enforce posi- 
tivity bounds [2(| on Sivers functions treating separately 
"valence" Sivers functions (u v and d v ) and "sea" Sivers 
functions (u and d). 



TABLE I. Best values of the free parameters for the Sivers 
function from fit to SIDIS data 0,0] on A$$* h ~' l " ) . 



X 2 /d.o.f. = 1.04 




= 0.05 ±0.05 


Ctd v 


= 0.76 ± 0.04 




= 0.78 ± 0.46 




= 2.09 ±0.25 


N Uv 


= 0.34 ±0.01 




= -1 ±0.01 




= fixed 


/3s e a 


= fixed 




= 0.003 ± 0.012 




= -0.15 ±0.01 


Ml 


= 0.45 ±0.10 (GeV/c) 2 







> 
3 



-0.05 - 



— l£ "D 0.05 

X 



13 



o 

0.01 


-0.01 



0.01 

IT3 o 

-0.01 




10" 



10" 1 



Fitting the pion data from both HERMES and COM- 
PASS we obtain a very good description of SIDIS data, 
with x 2 /d.o.f. = 1.04. The resulting set of parameters 
are presented in Table fl] together with the corresponding 
errors. As one can see, the biggest uncertainty is on pa- 
rameters P Uv and (3d v ■ This happens because SIDIS data 
covers a rather limited kinematic region in x < 0.3, as 
seen clearly in the HERMES plot Fig. [TJ Note that fu- 
ture measurements of JLab 12 [23j will explore the high- a; 
region in SIDIS which is very important for jet An as far 
as integration over x is performed in Eq. (J3]). 

In order to find the region of allowed values of f} Uv 
and Bj, we perform the scan procedure, also used in 
Ref. [24| to study the Collins effect. We produce a grid 
of values B Uv , Bd v G [0,4] in steps of 0.25 and for each 
pair of j3 Uv , Bd v perform a fit of SIDIS data. The re- 
sulting sets of parameters corresponding to 289 pairs of 
B Uv , Bd v give very good description of SIDIS data with 
X 2 /d.o.f G [1.04,1.08]; they are all almost statistically 
identical. Using these 289 sets of parameters we draw 
the shaded corridor in all the plots. It is important to 
realize that this corridor corresponds to almost the same 
description of SIDIS data. 

We present a comparison to the SIDIS data in Fig. [TJ 
which gives a very good description of HERMES tt + 
data. For ir°, ir~ asymmetries and zh and Ph± depen- 
dencies (and COMPASS data), the description is similar. 
In Fig. [2] we show the first fcj_-moment of the extracted 
quark Sivers functions versus x. 

We now assess whether the recently measured jet spin 
asymmetry from the AnDY experiment is compatible 
with the SIDIS Sivers asymmetry data; in other words, 
whether the jet asymmetry is consistent with our expec- 
tation on the process-dependence of the Sivers effect. To 
this end, we calculate the jet asymmetry An from Eq. (|4]) 
with our 289 equally-good sets of parameters. The result- 
ing shaded region for jet An as a function of Feynman 
xp = 2Pj z /yfs (Pjz the jet longitudinal momentum) is 
shown in Fig. [3J We note that the preliminary jet data 



FIG. 2. The first fcx-moment of the quark Sivers functions 
as a function of x, here /fj (a?) = —T q ^{x,x)/2M. Dashed 
lines correspond to positivity bounds, while the solid lines and 
the shaded region are the same as in Fig. Q] 
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FIG. 3. Description of AnDY preliminary data [l(| for inclu- 
sive jet production at forward rapidity (y) — 3.25 and center- 
of-mass energy y^s^OO GeV. Shaded region corresponds to 
the parameter scan. 

are inside the shaded region, which demonstrates that 
SIDIS Sivers data and jet An data are statistically com- 
patible with each other and that there exists a set of 
Sivers functions which can describe them simultaneously. 
Although we cannot claim that our analysis proves the 
process-dependence of the Sivers effect (the role of ISIs 
and FSIs) due to the large uncertainty and the very small 
size of the jet asymmetry data, at the very least such a 
process-dependence is not in disagreement with the ex- 
isting experimental data. Thus, we conclude that this 
is the first indication for the process-dependence of the 
Sivers effect. 

The very small size of the jet asymmetry is largely due 
to a cancellation between u and d quark Sivers functions, 
which have opposite signs. In order to carry out a more 
definitive test on the process-dependence of the Sivers 
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FIG. 4. Prediction of direct photon Am in pp collisions at 
rapidity y = 3.5 and center-of-mass energy ^=200 GeV. 



function, and simultaneously explore the consistency of 
the TMD and collinear twist-3 factorization formalisms, 
it is advantageous to study hadronic processes that have 
larger spin asymmetries. In this respect DY production is 
the ideal process to explore process dependence, provided 
that the effects of TMD evolution are completely under- 
stood; while direct photon production (though experi- 
mentally challenging [26j]) can be used to study the con- 
sistency of the factorization formalisms. First we make a 
prediction for the spin asymmetry An for direct photon 
production at RHIC kinematics, y — 3.5 and y/s = 200 
GeV, in Fig. |4] Since u and d quark Sivers functions are 
now weighted with their electric charge squared, which 
compensates the cancellation between them, we found 
that the direct photon An has much larger size ~ 5%, 
and it is negative 27J due to the nature of ISIs associated 
with the Sivers effect for direct photon production. 

Now, using the TMD factorization formalism [3(|, we 

sin(0-,— 4>s) 



UT 



compute the DY spin asymmetry An = — A 
at center-of-mass energy -/s = 500 GeV, for invariant 
mass 4 < Q < 8 GeV and transverse momentum < 
q± < 1 GeV. Due to the same reason, the asymmetry is 
larger ~ 8% and negative [25| (see Fig.[S]). At small and 
intermediate Xf region, the behavior is very similar to 
those in [3(j. The large xp behavior is different though 
it has a similar uncertainty band. 

In summary, we have analyzed the SSA for inclusive jet 
production in pp collisions collected by AnDY experiment 
and the Sivers asymmetry data from SIDIS experiments. 
We study the effect of ISIs and FSIs between the active 
parton and the remnant of the hadron on the process- 
dependence of the Sivers effect for both processes. After 
carefully taking into account all ISIs and FSIs in the in- 
clusive jet production, we find that the calculated jet 
spin asymmetry based on the Sivers function extracted 
from SIDIS experiments are consistent with the AnDY 
experimental data. Our result provides a first indica- 
tion for the process-dependence of the Sivers effect and 
further demonstrates consistency between the TMD and 



FIG. 5. Prediction for Sivers asymmetry An = — A^ 1 ^ 7 < ^ >s ' ) 
for DY lepton pair production at yfs = 500 GeV, for the 
invariant mass 4 < Q < 8 GeV and transverse momentum 
< q± < 1 GeV. 



collinear twist-3 factorization formalisms. However, due 
to the large uncertainty of the current preliminary data 
from AnDY and the small size of the jet spin asymme- 
try, our result cannot provide conclusive confirmation for 
process-dependence. Thus we also propose direct photon 
spin asymmetry along with DY measurements to test the 
process dependence of the Sivers effect. They are comple- 
mentary to each other, DY production is the ideal process 
to explore process dependence, provided that the effects 
of TMD evolution are completely understood; while di- 
rect photon production can be used to study the con- 
sistency of the factorization formalisms. We make pre- 
dictions for these asymmetries at RHIC kinematics, and 
encourage RHIC to undertake both measurements. 
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DE-AC05-06OR23177 (A.P.). 
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